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The reasons for the feasibility of the Microcanonical Fermionic Average (MFA) approach to lattice gauge theory 
with dynamical fermions are discussed. We then present a new exact algorithm, which is free from systematic 
errors and convergent even in the chiral limit. 



1. INTRODUCTION 

The dynamics of a gauge-fermion system on a 
d— dimensional lattice is described by its action 



S = -fiSpi(U) - Nf In det A(m, U), 



(1) 



where the fermionic degrees of freedom have been 
integrated out, and Nf is number of flavors tak- 
ing into account the species doubling. S p i is the 
plaquette energy [S p i = Y, p Re(U p ) / 'N c , for a 
compact theory, and S p i = — J2p®p/2, f° r non- 
compact QED]. 

The fermionic determinant det A (to, U) ac- 
counts for the effects of dynamical fermions; its 
numerical evaluation requires in general a huge 
CPU time due to nonlocality and dependence on 
different bare parameters. 

Therefore, one of the most challenging tasks in 
lattice simulations is the inclusion of dynamical 
fermions. Furthermore, most of the conventional 
algorithms can not be applied to the chiral limit 
(to = 0) because of the problem of convergence. 

In this paper, we discuss two efficient ap- 
proaches to simulating dynamical fermions which 
are convergent even in the chiral limit: the Mi- 
crocanonical Fermionic Average (MFA) method 
and a new exact algorithm. 
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2. FEASIBILITY 
METHOD 



OF 



THE 



MFA 



The main point of the MFA approach ffl to 
dynamical fermions is to compute the full effec- 
tive action as a function of the pure gauge energy 
E and bare parameters (to, Nf, and (3): 

S eff (E,m,N f ,p) = — In N(E) - ^—^lv/3E 



+ SF f (E,m,Nf), 
where 

N(E) = f [dU]6(S pl (U) 



d(d-l) 



VE) 



is the density of states at fixed E and 
SffJE, to, Nf) = — ln([det A (to, U)] N f ) 



(2) 
(3) 

(4) 



is the effective fermionic action computed as a 
microcanonical average over the probability dis- 
tribution 



S(S pl (U) - l^VE) 
N(E) 



(5) 



Then the partition function becomes a one- 
dimensional integral 



Z= / dEe- s °^ E ' m ' N f'P\ 



(6) 



Once the full effective action as a function of E 
and other bare parameters is known, the dynam- 
ics of the system can be investigated numerically 
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or analytically. There are several advantages in 
this approach: 

(a) S^ff does not depend on 0; 

(b) The dependence on Nj 1 is trivial; 

(c) can be simultaneously computed for any 
m (and particularly for m = 0) if the Lanczos 
algorithm is employed to evaluate det A. 

Now the most important question is how to 
compute S^ff defined by (4). Since the config- 
uration is sampled with the flat distribution (5), 
it may require high statistics to have the statisti- 
cal errors in det A under control. A more reliable 
way is to perform a cumulant expansion, which 
turns out to be an expansion on Nf 

-S[ ff (E,m,N f ) = AT/- (In det A) E 

-f-(Af L ) 2 ((lndetA- <lndetA) B ) 2 ) B + ..., (7) 

where the first term is some kind of mean field 
approximation, and the fluctuations correspond 
to the first corrections. 

As is well known, the mean field approximation 
is a good one if the interaction is long ranged or 
the dimensionality is large. In 0], we have given 
explicit physical examples (the local Ising model, 
the nonlocal Ising model, the pure gauge, mixed 
compact-noncompact model etc.) to show how 
the locality and dimensionality play a role. In the 
nonlocal Ising model, only the first contribution 
to the cumulant expansion exists. 

In the most interesting application of the MFA 
method, the gauge- fermion system, because of the 
nonlocality of det A, it is expected that the cumu- 
lant expansion (7) converges rapidly, which has 
indeed been observed in all the simulations P-[Io[| 
of lattice gauge theory with dynamical fermions 
(QED 2 , QED 3 , QCD 3 and QED±). For instance 
]2j, in noncompact QED4, the first cumulant is 
dominant, the second one is 2 % of the first one, 
and the third one is compatible with zero. Using 
the MFA method, the results for the Schwinger 
model @ and for QED A [| are in agreement 
with the exact calculable ones and consistent with 
those obtained with the Hybrid Monte Carlo al- 
gorithm (HMC) respectively. The simulation of 
QCD4 is in progress. 



3. EXACT ALGORITHM FOR THE 
CHIRAL LIMIT 

As mentioned previously, the feasibility of the 
MFA method depends on the convergence of the 
cumulant expansion. In practice, only the first 
few cumulants can be reliably calculated and a 
truncation of the expansion has to be made. Since 
a good control of systematic errors is difficult, 
especially when Nf is large, it seems interesting to 
us to develop some exact algorithm for obtaining 
quantitative results. 

An efficient exact algorithm for dynamical 
fermions at m ^ does exist, i.e. the HMC 
algorithm. In order to study the chiral proper- 
ties of lattice gauge theory without introduction 
of the external field rwi/jip or mass extrapolations, 
it is important to construct an exact algorithm 
which is convergent even in the m — » limit. To 
our knowledge, the Lanczos algorithm is the only 
efficient technique for evaluating det A which is 
convergent in this limit. The problem with this 
algorithm is that, if the canonical simulation is 
carried out by using the probability distribution 
e -(s g -A^indctA) for each Unk update thcn thc 

computer time is prohibitive; if many links are 
updated before testing for acceptance, then its 
rate is very low. 

Here we present a new exact algorithm, based 
on the following decomposition of the full action 

S = S g + S? f f + {-Nf In det A - S? ff ), (8) 

where S^f t is a local effective action chosen to im- 
prove the acceptance rate of the canonical simula- 
tion, and the partition function can be rewritten 
as 

Z = J [dU]e-( S * +S *lf^~(- N ? mdet A-Sf„). (9) 

This exact algorithm is implemented as follows. 

(a) Generate a new configuration U' with the 
probability distribution e~( S3+Se ff>, which can 
be done efficiently and is equivalent to a effec- 
tive pure gauge theory if we choose S^ft as that 
in Sect. 2. 

(b) Accept U' according to the probability 

Nf In det A(m,U')+S* ff (m,U') 

PA=min{l, eJV/MndetA(m ^ )+gf//(m ^ ) }, (10) 
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Table 1 



Acceptance rate of the exact algorithm 









sf„ = o 




S* ff obtained by MFA 










m = 


m = 0.02 


tyl = m = 


0.02 


N f 


= 2,f3 = 


0.3 


26% 


36% 


40% 


60% 


N f 


= 4,(3 = 


0.25 


14% 


15% 


26% 


30% 



which is just the Metropolis test. 

It can be shown that this approach is free from 
systematic errors and that detailed balance is sat- 
isfied. Furthermore this algorithm is convergent 
even in the chiral limit, since det A is evaluated 
by the Lanczos method. 

According to (10), the acceptance rate is higher 
if S^fj is a good local approximation to the 
fermionic action. Here we choose as that 

in the MFA method discussed in Sect. 2. 

We have tested these ideas in noncompact 
QED% using this new exact algorithm. To avoid 
huge changes in In det A, we change only part of 
the links in each iteration (a). The results for 
the 10 3 lattice (only 1000 links are systematically 
changed when U' is generated) are the following: 
(i.) For quenched (Nf = 0) and Nf = 2 theories, 
the data for the chiral condensate ('0V') at nonzero 
mass agree with those jllj obtained by the hybrid 
algorithm. 

(ii.) The results for Nf = 4 are slightly different 
from those in JHJ . This is not surprising because 
systematic errors increase as O(Nj) in the hybrid 
algorithm. 

(iii.) In the canonical simulations, we find that 
the choice of S^ff plays an important role in 
the improvement of the acceptance rate. Ta- 
ble |l| shows the results for the acceptance rate 
for different bare parameters. The data for direct 
canonical simulations (S^ff=0) are also included 
in this Table. As one sees, the use of the effective 
fermionic action obtained by the MFA method 
gives much better results for this quantity (it in- 
creases by about 40%). 

4. CONCLUSIONS 

To summarize, we have argued that one of the 
most important reasons for the feasibility of the 
MFA approach is the nonlocality of det A. Based 



on the MFA method, a new exact algorithm, 
which is convergent even at m = is presented 
and possible ways of increasing the acceptance 
rate are discussed. An interesting application of 
this exact algorithm, the investigation of the vac- 
uum chiral properties directly at m = by calcu- 
lating the probability distribution function of the 
chiral condensate, can be found in fl2fl . 
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